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CN ■ Motivated by the recent progress in cooling and trapping polar molecules, we present a simplified 

version of the rigorous contact pseudopotential for anisotropically-interacting polarized particles 

G\ ; [A. Derevianko, Phys. Rev. A 67, 033607 (2003); Phys. Rev. A 72, 039901(E) (2005)]. The 

simplifications are carried out for a practically important case of harmonically confined particles 
p ■ described by sufficiently smooth wavefunctions. The resulting contact pseudo-potential depends 

r^ , on the K-matrix of the underlying scattering process and is represented as a sum over pairs of 

C^ ■ 

^ I partial waves coupled by the collision. The contribution of each pair of the coupled waves {£ and 
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i' ) involves a tensor product of derivatives of orders I and i' . The asymmetry in the derivatives 
reflects the anisotropy of the original interaction potential: there is a preferential appearance of 
the derivatives along the polarizing field. 
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I. INTRODUCTION 

The intricacies of the many-body problem are rooted in inter-particle interactions that 
lead to non-separable Hamiltonians. Introducing pseudo-potentials, i.e., effective interac- 
tions that are simpler than the original interactions, makes the problem more tractable. In 
particular, due to their compelling efficiency, the contact (zero-range) pseudo-potentials p| 
have enjoyed a remarkable success in a number of applications, the physics of ultracold 
atomic gases being one of the most recent examples. In particular, properties of Bose- 
Einstein condensates may be well understood just in terms of a contact potential (Fermi [2]), 
with its strength determined by a single scattering parameter - scattering length |3|, |j]. 

Due to their importance in ultracold atomic physics, there has been a recent resurgence 
of_ interest in the pseudo-potential approximation, see, e.g., Refs. |5|, Ig, |7|, |8|, |9|, 
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12|. All of the enumerated works focused on particles that predominantly interact via 
spherically-symmetric forces. Here, motivated by recent experimental developments, we 
focus on anisotropically interacting atoms and molecules. Indeed, much of the present 
experimental efforts on cooling and trapping of atomic particles turned to "unconventional" 
(non-alkali) species: molecules [l3| and atoms with complicated open-shell structure [ij, 
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16| . Most of these species possess internal multipolar moments, such as electric and 
magnetic dipoles, quadrupole moments, etc. As a result, polarized atoms and molecules 
interact via anisotropic forces: the interactions depend on orientation of the inter-particle 
axis with respect to the orientation of the internal moments. 

Previously, I derived a rigorous contact pseudopotential for anisotropically interacting 
species [l7|. While being general, the rigorous prescription is hardly necessary for smoothly 
behaving wavefunctions (mathematically, we require the existence of the Fourier transform 
(i.e., the momentum representation) of the wavefunction and all of its derivatives). More- 
over, most of the cold-atom experiments are carried out in a presence of harmonic trapping 
potentials. These two observations facilitate simplifications of the rigorous pseudopotential. 
In particular, we show that integro-differential operators entering the rigorous pseudopoten- 
tial may be replaced by a mathematically simpler combination of derivatives. 

One may distinguish between two classes of applications of pseudopotentials: (i) deter- 
mining exact few-body solutions and (ii) usage of pseudopotentials within a many-body 
mean-field framework. In the former case one needs to carefully deal with regularization 



issues and the rigorous pseudo-potential is required. In the later, many-body, case, the 
functional space is restricted and the regularization is not needed. The results of this paper 
are geared towards the many-body applications of the pseudopotential method. 

This paper is organized as follows. First, In Section [III we formalize the problem and 



recapitulate results for the rigorous pseudopotential of Ref. [l7|. In Section [XTTl we consider 
simplifying conditions of the collision process in a harmonic trap. Before proceeding to gen- 
eral derivation in Section|Vl in Section[IV]we illustrate the involved techniques by simplifying 
a truncated dipolar pseudopotential and the contact Fermi interaction. In Sections I VI Al and 
IVIBl we consider two limiting cases of the isotropic interactions and anisotropic interactions 
in the Born approximation. Finally, the summary is given in Section IVIII 

II. PROBLEM SETUP 



While the formal setup of the problem is identical to that of Ref. [17|], below I would 
like to recapitulate several points as they relate to more recent developments in ultracold 
physics with polarized heteronuclear molecules. 

In a typical cold-molecule setup, a dilute gas of polar molecules forms a cloud in an exter- 
nal harmonic trapping potential. Orientation of molecular dipoles D is fixed by application 
of a polarizing electric field (otherwise, molecular rotations would average dipole moments 
to zero). At large inter-molecular separations, the molecular interactions have the dipolar 
character, 

V{j)^^(?>cosH-l). (1) 
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Here 6 is the angle between the inter-molecular (collision) axis r and the polarizing field (it 
is assumed to be directed along the z-axis). Clearly the interaction is anisotropic: if the 
collision axis is aligned with the external field, the dipoles attract each other, and if r is 
perpendicular to the field, the interaction is repulsive. 

Dipolar interactions are also important for the BEC of highly-magnetic chromium 



atoms 18|. Again the species are polarized by applying an external (magnetic) field. The 
atoms are trapped in the lowest-energy Zeeman sublevel; the depolarizing transitions to the 
upper-energy levels are forbidden energetically. 

While most of our discussion focuses on the dipolar interactions, the final result will be 
also valid for a general case of axially-symmetric interactions. We will parameterized such 



interactions as 

\/(r) = VL(r)PL(cos ^) , L is even , (2) 

where Pl is the Legendre polynomial. Due to the axial symmetry, L is an even integer. In 
particular, for isotropic interactions, L = 0. For dipolar interaction ([T]), L = 2, V2(r) = 
-2D^/r\ 

Constructing the pseudopotential in the Huang- Yang approach [l9| , essentially boils down 
to finding an effective contact interaction which, at sufficiently large separations, reproduces 
the exact solutions of the Schrodinger equation with the full original potential. To this end 
one requires knowing the exact solutions of the scattering problem with the full potential 
in the asymptotic region. The relevant radius of the asymptotic region may be estimated 
by comparing centrifugal term in the corresponding radial equations with the multipolar 
contribution. Parameterizing Vi(r) = C^/r^, the radius is 

Rk- [ TjCk 1 • (3) 

While the pseudopotential would reproduce the solutions for r ^ ^fc, it is worth empha- 
sizing that the scattering solution depends on the entire interaction potential: as molecules 
approach each other, the electronic clouds start to overlap, and the interactions substantially 
depart from the dipolar form ([1]). In addition, the anisotropic interactions couple various 
partial waves: for example, the dipolar interaction nominally couples partial wave i with 
partial waves i ±2. In other words, one has to solve a multi-channel scattering problem. 

The exact linear combinations of the regular and irregular free-particle solutions in the 
asymptotic region are naturally expressed in terms of the K-matrix (reactance matrix), 
^im 5 which characterizes couplings between {£m) and {£'m') partial waves. To formalize 
the concept of the K-matrix, we consider scattering wavefunction $ (r) for a pair of particles 
colliding with the relative linear momentum of hk. In the asymptotic region it may be 
expanded using free-particle solutions 

$ (r) ^ ^ {aimji (kr) - Pi^rii (kr)) Yi^ (r) , (4) 

Im 

where ji [kr) and rii {kr) are the spherical Bessel and Neumann functions respectively and 
aim and jSim are integration constants. These constants are related requiring the solution $ 



of the Schrodinger equation with full potential to be regular at r = 0, 

Am = ^ /C;^' avm' ■ (5) 
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This relation defines the K-matrix. 



The anisotropic pseudo-potential [l7| reads 



V'-nr) = 4 E ^^^«". (f) % ("'"»'*) . (6) 
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where 5 (r) is the Dirac 5-function and the projection operator is defined as 

. f (j \ ^''+^ f 

The matrix element of the pseudopotential ([6]) between two plane waves (r|k) = 
(27r) e^^'^ is given by 

ll';mm' ' 

In the following we assume that the global many-body properties of the quantum gas can 
be described by well-behaved wavefunctions. We will base our derivation on the momentum- 
space kernel Eq. ([8]). This leads to the requirement that the momentum representation 
of the wavefunction and of all of its derivatives must exist. This constraint follows from 
the fact that one needs to interchange differentiation and integration operators in Eq.© 
while arriving to Eq.Q for a momentum-space wavefunction. This constraint, for example, 
discards wavefunctions diverging as 1/r at the origin. 

Let me explicitly address a somewhat confusing point of dipolar interactions being of 
"long-range" character. Can the contact pseudopotential approximate an action of the 
"long-range" dipolar interaction? The use of the contact pseudopotential is mathematically 
justified because there is a radius ^3, Eq.( [3]), where the centrifugal forces overtake the 
dipolar interaction. Beyond ^3 the wavefunction no longer "feels" the dipolar interaction, 
as reflected by Eq.( Hj). In this sense, the dipolar interactions are not "long-range" but 
rather "long(ER)-range" . The true anisotropic "long-range" potentials would have \jl^ 
tail, where k < 2. Notice that in this discussion, I have used the fact that the anisotropic 
potentials do not directly contribute to diagonal part of the i = radial equation void of 
the centrifugal contribution. 



III. SIMPLIFICATIONS FOR HARMONICALLY-CONFINED PARTICLES 

Now, under simplifying assumption of harmonic trapping, I transform the momentum- 
space expression ([H]) back into the coordinate space. We write for a matrix element of the 
pseudopotential (cf. Ref. [20| for Rydberg atoms) 

(^iKpslV^) = (27r)"^ I dkdk'drdr' X 

t/j* (r) e-'^-^v (k, k') tp* (r') e^^'-"' . (9) 

Below we simplify this expression by noticing that only certain values of |k| and |k'| con- 
tribute to this integral. 

Experimentally, the collisions occur in the presence of a harmonic trapping potential 

u (r) = ^M {uy + uy + uy) , (lo) 

where Ui are the trapping frequencies. For two harmonically-confined particles the center- 
of-mass and relative motions decouple and the Hamiltonian for the relative motion reads 

where V [ru) is the full interaction potential between the particles and fi = M/2 is the 



ue equation Hrip (ri2) = 
presented a numerical 



reduces mass. In the stationary problem, we solve the eigenva^ 
Erip {ri2), Er being the energy of the relative motion. Ref. p, u\ 
comparison of solutions of this equation with the full molecular potential and its pseudopo- 
tential representation (for isotropic s-wave scattering). Their main conclusion is that the 
pseudopotential description holds as long as the range of the potential is much smaller than 
the oscillator length Li = {Ti/MuJiY/'^. Ref. [7| also observed that for velocity-dependent 



potentials the relevant collision momentum in Vps is ?i fc^/(2yu) = E^. Indeed, the collision 
process occurs at |ri2| much smaller than the harmonic length. In this region [/ (|ri2|) ~ 
and the kinetic energy is E^. 






Recently, there was a study [211] of the validity of the anisotropic pseudopotential 
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approach 



for dipolar interactions [17| in a harmonic trap. Similarly to the isotropic casep, l7|, these 
authors found that the pseudopotential description remains accurate as long as the dipolar 
length ^3 is smaller than the characteristic length of the trapping potential. Their numerical 
experiments also support the discussed choice of the collision momentum kc- 



Based on the above comparison with the exactly-solvable two-body case, we require that 
the multipolar lengths, Eq. ([3]), are much smaller than the harmonic lengths, Rl+i <^ niin Li, 
i.e., the pseudo-potential treatment is valid in the region where the wavefunction remains 



ude 



2l|. 



essentially flat. We notice in passing, that the more rigorous consideration should inc 
criteria related to the scattering parameters of the problem, such as the scattering length 
Under this condition, the interparticle interactions have a negligible effect on the energies 
of the harmonic motion, as the characteristic multipolar interaction energy is much smaller 
than the excitation energy of the harmonic potential and we identify E^ ~ h{nxUJx + UyUy + 
rizUJz + 3/2), rii = 0, 1,2, . . .. For the ground state of the harmonic oscillator the above 
discussion leads to the value of the collision momentum of 

kl^3/2Muj/n, (11) 

with u! = X]i^«/3 being the average of the three trapping frequencies. The generalization 
to the excited states is straightforward. 

We also notice that for a many-particle case, due to mean-field effects, the effective 
potential felt by the particles may deviate from the external harmonic trapping potential, 
Eq. flTOl) . In this case the mean- field wavefunction needs to be fitted to a corresponding 
harmonic oscillator wavefunction and the kc has to be consistently redefined based on the 
fitting. Indeed, most of the variational studies of BECs employ the Gaussian ansatz for 



the trial wavefunction, see for example, Ref. 22]. In this sense, the choice of k^ would be 



self-consistent in such variational studies. Otherwise, the suggested fitting is an approximate 
prescription. 

Returning to the evaluation of the integral ([9]), we see that the relevant contributions are 
accumulated at |k| = |k'| = k^. In other words, we consider only the "on-shell" collisions, 
with the radius of the shell determined by the energy of the harmonic motion. This situa- 
tion is similar to the interaction of Rydberg electrons (of well-defined kinetic energy) with 



atomic particles, considered by Omont 20|]. The present derivation may be considered as a 
generalization of his approach to anisotropically interacting particles. We will recover the 
Omont's result as a limiting case of isotropic interactions in Section fVI A[ 



IV. TRUNCATED DIPOLAR PSEUDOPOTENTIAL 



Before presenting the general analysis, I will demonstrate the procedure and discuss some 



subtleties for a simpler case of the truncated dipolar pseudopotential [17|]. In this model, 
the entire summation over the coupled partial waves in the momentum-space kernel (jH]) is 
limited to the s — s and s — d couplings, 

1 h"^ 

■V (k, k') ~ ^^ i^ss - asd T (k, k')) , 

with 

T (k, k') = VS {P2 (cos Qk) + (fc/A;')' P2 {cos 9k'] 

where dk and 6k' are angles between the polarizing field and k and k'. Parameter ass is the 
conventional scattering length 

ttss = ~'^00v"'//"' ' 

and ttsd is the so-called "off-diagonal" scattering length, characterizing the strength of 
anisotropic coupling between the s and d waves 

ttsd = -K?^{k)/k. 

As A; — i> 0, both expressions approach finite values. 

Let us consider the isotropic (a^s) term first. We expect to recover the original Fermi 
contact interaction, Anh? /M ass ^ (r), as the result of the analysis. 

^ ^ rfkrfkV^'^'t; (k\ k) e-^'^" = \ \ ^Qss /rfk [ dk'e'^'''' le"''^-'" = 



{27Tf J ' ' ' (27r)'27r2M 



-^—ass(27i?S(r)6(r') 
3 



where we used the identity / d^k exp[ik • r] = (27r) 6 (r). Then 

I h^ .. ,3 f f I n"" 



I^P.|^) = ^M<'ss {2nr / 6{r)6 (r') r (r) ^ (r) drdr' = ^^a.s (27r)^ r (0) ^ (0) 



From here we read 



Vps{r) = An—assS{r) , (12) 



which is, indeed, the Fermi contact interaction. 

Now we proceed to simplifying the dipolar (asd) term. We focus on the following integral 

rfk / dk'e"^'-"' P2{cos9k)e-'^-'' . 
8 



The resulting operator form of the pseudopotential depends on the choice of the representa- 



tion of the Legendre polynomial P2 (cos 6k) '■ 
as 

P2icosek] 



cos^ Ok- We may express the polynomial 

3fc| 
' 2A;2' 



or, equivalent ly, as 



Po (cos I 



11 'Mc'^ 



^5(r)5(r)-^^5(r)^5(r) 



2 p v-x . -y ' -zi . 2 A;2 • 
While computing the integral, we evaluate the |k|-dependent functions at the "on-shell" 
value of |k| = kf. and pull them out of the integration. We also use J/c^exp(ik • r) rfk = 
(l) ^ /exp(ik ■ r) (ik = {2n) ^^5 (r). With the first combination we obtain 

I = 

and with the second combination 

Both combinations are equivalent for harmonic-oscillator wavefunctions. The equivalence 
depends on the definition of the collision momentum k^, it has to be chosen consistently. 
In the first case we arrive at the pseudopotential 



Ks (ri2) 



4:71 ^6 (ru) X 






< ttss + V^ttsd 
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^ + 2^2 


f)2 ' f)2 



(13) 



and in the second case to the equivalent pseudopotential 

Ks (ri2) ~ 47r— 5 (ris) x 



M 



(14) 



+ ^asd-rr; 



kl 
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In these expressions, the operator ^ acts on the bra and ^ operates on the ket. In both cases 
the dipolar contribution to Vps breaks into two parts: isotropic and anisotropic terms. The 
anisotropy is imposed by the polarizing field (the derivatives in the anisotropic contributions 
are taken along the field). The long-range character of the dipolar interactions is manifested 
through k\ which characterizes the entire trapping potential. 



Being equivalent, which representation, Eq. (fT3!) or Eq. (fT^ . is more accurate in approx- 
imate calculations? A similar question appears in choosing the gauge-dependent operators 
(e.g., velocity-gauge and length-gauge of the dipole operator) when carrying out compu- 
tations involving radiation fields in atomic physics or quantum electro-dynamics. When 
carrying out approximate calculations the results would differ. The advantage of repre- 
sentation (fT3l) is that the isotropic part of the dipolar operator participates on an equal 
footing with the Fermi contribution: the isotropic dipolar contribution merely renormalizes 
the traditional spherically-symmetric pseudopotential; this may be beneficial in interpreting 
the results of the calculations. The advantage of the second representation ( TT^ is that it 
is the zeroth component of the rank 2 irreducible spherical tensor (we will discuss this in 
detail in the following Section). As such its expectation value manifestly vanishes for any 
spherically-symmetric wavefunction, while the expectation value of flT3|) vanishes only for 
the spherically-symmetric states of the harmonic oscillator. 

Since the properties of dipolar quantum gasses are usually studied as a function of the 



aspect ratio of the trapping potential, see e.g., 23|, l2J], the tensorial representation seems 
to have a computational advantage. Another advantage is that the resulting expressions for 
higher partial waves are more concise. Due to these considerations, below I will treat the 
general case using techniques of irreducible tensor operators. 

V. GENERAL CASE 

In this Section we carry out simplifications of the rigorous pseudopotential ([6]) for an 
arbitrary potential and all partial waves. The derivation essentially follows the steps for the 
truncated dipolar pseudopotential of the preceding Section. I will use several formulae from 
the quantum theory of angular momentum. The notation and phase convention follow that 



of the standard compilation [25|]. In particular, normalized spherical harmonics (C-tensors) 
are defined as C^(f) = \/^fi ^m(r)- The tensor product of irreducible tensor operators 
(such as C-tensors) is expressed in terms of their spherical components as 

1^00 ^ 50.)}^^^ = Y: Ci^.,^,J^^ , (15) 

mim,2 

where C'j^mii2m2 ^^^ ^^^ Clebsh-Gordan coefficients. Beyond streamlined notation, there are 
several advantages to operating in terms of the tensor products, e.g., angular selection rules 

10 



and a suitability for applying the Wigner-Eckart theorem. We will also make use of the 
following representation of the C-tensor 

p,q,r -f^ ^ \ / \ / 

where p, q, and r are non-negative integers satisfying p + q + r = I and p — q = m. 

In addition we will use general properties of the K-matrix }C\^ . First of all, the polarizing 
field imposes the axial symmetry on the collision process. This symmetry leads to the 
conservation of the projection of the angular momentum on the direction of the polarizing 
field (z-axis) and to the selection rule m = m'. Parity conservation requires that only the 
partial waves of the same parity are coupled: for example, the p-waves are uncoupled from 
the s and d waves. Formally, 

We also relate /C^ 'Zj^ = /C^ '^J^ . Finally, the time-reversal symmetry leads to ICf^' = K}j?^, . 
Again it is worth reminding that we limit our consideration to axially-symmetric components 
of the full potential ([2]), where L is even, and the enumerated properties can be derived by 



a direct examination of the radial scattering equations (see, e.g., appendix of Ref. 17l|). 

We may rewrite the momentum-space kernel, Eq. ([8]) as an expansion over tensor products 
of normalized spherical harmonics 



% (k, k') = 7^ E E ^^ (^^ ' ^^) ^''" {C' (k) ® C^' (k') },,,^„ ^ (17) 

J even IV 

with the coupling coefficient 

Tj {W; K) = -An^ [{21 + 1)(2/' + 1)]^/^ J] (-1)™ C£^=!^^^^^ • (18) 



In the sum, the compound angular momentum J runs over all even integers in the range 
\l — l'\ < J < l + l'. Since both / and /' are simultaneously even or odd (parity selection rule), 
the sum is limited to even values of J. In addition, following the arguments of Section IIIII 
we have evaluated the kernel at the on-shell value of |k| = |k'| = kc- Notice that only the 
M = component of the compound tensor participates in the expansion - this is due to the 
axial symmetry of the collision process. 

Now we turn to evaluating the integral over momenta in Eq.([9]). We expand the C-tensors 
according to Eq. (IT6|) and use the properties of the Dirac (5-function as in Section [IVl On 
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recombining the resulting expansion we arrive at 



or in the operator form 



J even W "" I \ "^ / V "" / J J,0 



J even U> I V V \ "" / ) .J,0 

Here and below the application of the (5-function is understood in the sense of the limit in 
Eq. ( TT9l) . The C-tensors C^ \y) ^^^ differential operators of order Z, defined as 

The spherical components of the nabla operator V are expressed in terms of its Cartesian 
components 

\f2 \dx dy) 

oz 
Explicit representation of several differential operators C^ ( ^ ) is given in Table [B 



(^, 


m) 


kc^ 




(0,0) 


1 




(1, 


0) 


dz 




(1, 


±1) 


1 

V2 


(=F9a; + i9j^) 


(2, 


0) 


H 


-A + 35f ) 


(2, 


±1) 


^\ 


/|(9^±i9j,)a^ 


(2, 


±2) 


1 

2V 


/|(9,±i5,)2 



TABLE I: Explicit form of operators C^ f ^ J for several values of £ computed using Eq. ([2T]) . We 
use a shorthand notation for partial derivatives, e.g., d1 = ^. 
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We also mention an alternative representation of the C-tensors, which may be more useful 
in some applications 



^"0^(^^^<-»^*^''*^'--*^} 



Im 



By limiting the summation to pairs of partial waves s — s (resulting J = 0,171 = 0), s — d 
and d — s (resulting J = 2, m = 0), we recover the result, Eq. (ITll) . for the truncated dipolar 
pseudopotential. 

We would like to emphasize that so far we have not specialized our consideration to a 
particular multipole contribution. In fact the original interaction potential may contain 
a sum over several multipolar contributions. In the following sections we specialize our 
discussion to two limiting cases: isotropic interactions and weak anisotropic interaction, 
where the Born approximation is valid. In the former case we recover literature results and 
in the latter (Born) case we show that only J = L component of the compound tensor 
remains in the general expansion (l20l) . 

VI. SPECIALIZED CASES 

A. Isotropic interactions 

In this Section we simplify the general pseudo-potential expression, Eq. (120|) . for the case 
of isotropic interactions. In this case the various partial waves decouple and the reactance 
matrix reads 

/C!!™) = 6w 6mm' tan 6, , 



iso 



Si being the conventional phase shifts. This relation can be easily derived by comparing 
Eq. (jl]) with the textbook expressions for isotropic scattering. The expression for the coupling 
coefficient T, Eq. (fT8l) . simplifies to 

M kc 

where we used ^^ (—I)'" C'/m-i -m — ^Jfi i~^) a/(2Z + 1). For isotropic interactions we deal 
with only J = (scalar) component of the compound tensor. This reflects the rotational 
properties (scalar) of the full isotropic potential. Further, we invoke the theorem for addition 
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of spherical harmonics 

I V '^c / 

and arrive at 

h v^/ , N / tan (5; A ^ , X / V V \ , , 

(Vi.),„ = 2.- ^(2/ + 1) (-^j S (r) P, ^- . -j . (22) 

Again the apphcation of the (5-function is understood in the sense of taking the hniit in 
Eq. fimi) . This is a resuh previously derived by Omont [20| in the context of interaction of 
Rydberg electrons with atoms (in that case, the reduced mass fi ~ rrie; in our case /i = M/2). 
The / = (s-wave) contribution in this sum is the Fermi pseudopotential. The / = 1 and 
I = 2 {p- and d-wave) terms were re-derived recently in Ref. 9|. 

B. Anisotropic interactions in the Born approximation 

The derived pseudopotential is based on a microscopic description of the collision process. 
The relevant K-matrix is to be determined from the full solution of the Schrodinger equation 
with the original potential. If the perturbation imposed by the interaction potential on the 
free motion of the particles is sufficiently weak, one may employ the Born approximation. 
Of course, the exact K-matrix is not guaranteed to be equal to its Born value. The Born 
approximation would break, for example, near resonances, found in dipolar scattering in 



Ref. 26|. Even without resonances, according to Ref. 27|, the Born approximation for 
dipolar scattering is valid only if R^/R^^^ <^ 1, where -Rvdw ~ lOObohr is the van der 
Waals length [R^ based on the van der Waals dispersion coefficient C^). This inequality 
easily breaks for polar molecules, where ^3 ~ 10^ bohr. 

Nevertheless, it is instructive to consider the pseudopotential in the Born regime. To 
the lowest order, the K-matrix in the Born approximation for anisotropic interactions was 
derived in Ref. [l7|. It reads 

feir'). =-^l\vHlm\C^o\l'm')^ (23) 

V / Born n 

POO 

4^) = / x''dxVLix)ji (kx) J, (kx) . (24) 



Here {lm\C^\l'm') = J dnYj*^{n) Cq {^)Yi'm' {^)- Invoking the Wigner-Eckart theorem 
for this matrix element of the C-tensor, we may introduce the reduced matrix element of 
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the K-matrix, (/||/C(^)||/'), via 

^i:r' = (-ir"| ^ ^ ^' |(^II^^^^II0, (25) 

-m r/i 

(/||/c(^)||r) = -^4f(/||c(^)|ir), (26) 

with the reduced matrix element of the C-tensor being conventionally expressed as 

(/||C(^)||/') = (-1)' y/{2l + 1) {21' + 1) I ^ ^ ^' I . (27) 



By substituting Eq. fl25|) into Eq. flTSl) . summing over magnetic quantum numbers, and us- 
ing the orthogonality relation for the Clebsh-Gordan coefficients, we obtain the simplified 
expression for the coupling coefficient 

■^ ^ ' ^ M V 2L + 1 J K 

Accordingly, only the J = L compound tensor participates in the pseudopotential 



P^^Born.L M ^^ ^ \ 2L + 1 J k, 



W ^ ' ' ' '^ 



For the specific case of dipolar interactions, L = 2, selection rules require I' = I or 
/' = / ± 2. The reduced matrix elements of the K-matrix are proportional to the wave- 

n 

vector (c.f. introduction of generalized scattering lengths for dipolar scattering in Ref. |17]). 
Explicitly, 



(/||i^(^=2)||/')B„,,,DD = -D^^k {l\\C'\\l') X <; 



[/(/ + !)]"' /' = / 

[3 (/ + !)(/ + 2)]"' I' = 1 + 2 . (30) 
[3/(/-l)]"^ I' = 1-2 



By limiting the summation to / = 0, /' = 2 and I = 2,1' = couplings in Eq. fl29p . we recover 
the truncated dipolar pseudopotential, Eq. flT^ . where we use the off-diagonal scattering 
length in the Born approximation [17 1 

(Meo™--^^^. (31) 

15 



Notice that (ctsd)Born is ^^ ^^^ same order of magnitude as the radius of the asymptotic region 
for the pseudopotential, ^3, Eq. ([3]). 

Due to the restriction to the J = L,M = component of the compound tensor, the 
transformational properties of the pseudopotential under rotations of the coordinate frame 
are the same as of the original V^-component of the potential. While we have proven this 
property in the Born approximation, it seems not to hold in the non-perturbative regime. 
For example, in the case of dipolar interactions, radial equations for the m = component 
of the scattering wavefunction involve couplings oi s — d, d — g, g — h, ... partial waves, 
while for the m = ±2 components, the s — d coupling is missing. A resonance in the s-wave 
channel would drastically affect the K-matrix for the m = channels, while the m = ±2 
components are not affected by this resonance. In this case, the Wigner-Eckart theorem for 
the K-matrix, Eq. fl25l) . does not hold and we no longer arrive at the restriction J = L in 
the pseudopotential. 

Finally, it is worth pointing out the connection of the derived Born formulae to the original 
potential, Eq.(l2]). This also provides another test of self-consistency of the derivation. To 
this end, we transform the original potential to the momentum-representation 



v{k',k) = ^ j dre-'^'-WLir) PL{f)e 



By expanding the plane waves into the spherical waves, evaluating angular integrals, and 
setting the |k| = |k'| = fee (on-shell values) we arrive at the pseudopotential kernel, Eq.®, 
with the elements of the K-matrix evaluated in the Born approximation. This means that 
following the derivation of this paper we would necessarily arrive at the pseudopotential in 
the Born-approximation, Eq. (!29l) . even discarding the microscopic description of the colli- 
sion physics. This is another test of self-consistency of the present approach in the Born 
approximation. We emphasize again that a description starting from Eq. ([2]) would miss 
deviations of the scattering matrix from it's Born value. 

VII. SUMMARY 

To reiterate, here we derived a contact pseudo-potential for anisotropically interacting 
particles trapped in a harmonic potential. Such interactions include, for example, dipole- 
dipole interactions between polarized molecules and atoms and subsume the conventional 
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isotropic interactions. Due to the well-defined collision energy and assumption of sufficiently 



smooth wave-functions, the previously derived [17|] rigorous, but complicated, pseudopoten- 
tial is simplified. An additional assumption is that the characteristic radius of the potential, 
Rk is much smaller than the harmonic length of the trapping potential. All of these as- 
sumptions are natural in the contexts of studies of many-body physics, e.g., BECs of dipolar 
gases. 

The simplified contact pseudopotential, Eq. (120|) . involves a sum over partial derivatives. 
The anisotropy of the interaction is reflected in the preferential appearance of derivatives 
along the polarizing field. Compared to the rigorous pseudo-potential of Ref. [l7|], the 
presented simplified pseudo-potential may prove to be more practical in applications to 
trapped quantum gases of polarized atoms and molecules. The advantage of the pseudo- 
potential description is that, being based on the microscopic quantum description of the 
collision process, it remains valid even in the non-perturbative regime, where the Born 
approximation breaks down. 
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